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I review all the new features of the HERWIG event generator which are relevant 
to Linear Collider (LC) physics starting from version 6.1. 

1 Introduction 

HERWIG is a general-purpose Monte Carlo (MC) event generator for high- 
energy processes ^ , providing a full simulation of hard lepton-lepton, lepton- 
hadron and hadron-hadron scattering and soft hadron-hadron collisions in a 
single package, comprising the following features. 

1. Initial- and final-state QCD jet evolution with soft gluon interference taken 
into account via angular ordering. 

2. Colour coherence of (initial and final) partons in all subprocesses, including 
the production and decay of Supersymmetric (SUSY) particles. 

3. Matrix Element (ME) corrections to the parton shower (PS) algorithm. 

4. Spin correlations in the decay of heavy fermions. 

5. Lepton beam polarisation in selected processes. 

6. Azimuthal correlations within and between jets due to gluon interference 
and polarisation. 

7. A cluster model for jet hadronisation (with optional colour rearrangement) 
based on non-perturbative gluon splitting and a similar model for soft and 
underlying hadronic events. 

The HERWIG source codes and related information can be found in ^ . 



2 Linear Collider physics with HERWIG 

Of particular relevance to the programme of future electron-positron LCs are 
items 3.-5. We will describe these in the following, after mentioning the Min- 
imal Supersymmetric Standard Model (MSSM) processes currently available 
in lepton-antilepton annihilations. 
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2.1 MSSM processes in £'^£ collisions (i = e,iJ.) 

Starting from version 6.1, HERWIG includes the production and decay of 
(s)particlcs, as given by the MSSM. The conventions used and a detailed de- 
scription of the implementation can be found in ^ . 

HERWIG does not contain any built-in models for SUSY-breaking scenar- 
ios. In all cases the general MSSM particle spectrum and decay tables must 
be provided just like those for any other object, with the following caveats: (i) 
SUSY particles do not radiate (which is reasonable if their decay lifetimes are 
much shorter than the QCD confinement scale); (ii) CP- violating SUSY phases 
are not included. A package, ISAWIG, has been created to work with ISAJET 
^ to produce a file containing the SUSY particle masses, lifetimes, couplings 
and mixing parameters. This package takes the outputs of the ISAJET MSSM 
programs and produces a data file in a format that can be read into HERWIG 
for the subsequent process generation. The user can produce her/his own file 
provided that the correct format is used. To this end, we invite the consulta- 
tion of the ISAWIG webpage^, where some example of input files can be found 
(the SUSY benchmark points recommended in Ref. ^ are also available) . 



Table 1. The MSSM processes via i f. anniliilatious implemented in HERWIG. 



IPROC 


MSSM processes 


700-99 


R-parity conserving SUSY processes 


700 


£'^£~ — > 2-sparticle processes (sum of 710-760) 


710 


£'^£~ —> neutralino pairs (all ncutralinos) 


706+4IN1+IN2 


£^£~ X1N1X1N2 (INI , 2=noutralino mass eigcnstatc) 


730 


£^£^ — > chargino pairs (all charginos) 


728-I-2IC1-I-IC2 


£'^£~ — > X1C1X1C2 (ICl,2=chargino mass eigenstate) 


740 


£'^£~ —^ slepton pairs (all flavours) 


736+5IL 


£+£- ^ Il,rII^r (IL = 1, 2, 3 for 1= e, A, f) 


737-F5IL 


£+£- £l£1 (IL as above) 


738-h5IL 


£+£- £l£*r (IL as above) 


739+5IL 


t£~ £r£*ii (IL as above) 


740-F5IL 


£+£- VLvl (IL = 1,2,3 for z7e,^V,?T) 


760 


— » squark pairs (all flavours) 


757-I-4IQ 


£"^£- -> qL.RolM (iq = 1...6 for q = d...t) 


758-F4IQ 


£'^£~ —> qLqt (IQ as above) 


759-F4IQ 


£'^£~ —> qLqR (IQ as above) 


760-F4IQ 


£'^£~ — > qrOr (IQ as above) 



In addition to the decay modes implemented in the ISAJET package, 
ISAWIG also includes the calculation of all 2-body squark/slepton and 3-body 
gaugino/gluino R-parity violating (RPV) decay modes (alas, RPV lepton- 
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gaugino and slepton-Higgs mixing is not considered). Moreover, the emulation 

of RPV processes is also a feature of the production stage. Tab. 1 illustrates 
all MSSM modes initiated by scattering that are available at present 

(version 6.4). Here, IPROC is the input label selecting the hard process. 



Table 1. Continues. 



800-99 


R-parity violating SUSY processes 




800 


Singl 


c sparticle production, sum of 810-840 


810 






(all neutralinos) 




810+IN 


£+£- 




, (lN=neutralino mass state) 


820 


e+r 


-X-e+ 


(all charginos) 




820+IC 




^ Xice; 


'", (IC=chargino mass state) 


830 






and £+£- -> efW' 




840 


i+r 




H°/A° and £+£' 


£+H- 


850 


i+i- 






860 


Sum 


of 870 and 880 




870 


i+r 




via LLE only 




867+3IL1+IL2 


i+r 


^j^Ll^rL2 (ILI. 2=1,2,3 for 


e,H,r) 


880 


i+r 


— > dd, via LLE and LQD 




877-|-3IQl+iq2 


i+r 


rflLlf' 


IL2 (IQ1,2=1.2,3 for 


d, s, b) 



Table 1. Continues. 



910-975 


Higgs processes 


910 


£+£- ueuth° +£+£-h° 


920 


£+l~ ueueH° +£+£-H° 


960 


£+£- ^ 


970 


£+£- Z°H° 


955 


£+r H+H- 


965 


£+£- ^ ^"/i" 


975 


£+£- A°H° 



2.2 ME corrections to the standard PS 

The study of top quark production and decay will be one of the main areas 
of research activity at future LCs. Hence, it is of paramount importance the 
availability of a MC event generator describing the dynamics of the e+e" 
tt bbW~^W~ process {mt = 175 GeV) with great accuracy. As the top 
quark is a coloured particle, of particular concern are higher order effects from 
QCD. Whereas those induced by the exchange of virtual gluons do not change 
the lowest order kinematics, real radiation of the latter does so. This is of 
relevance when it comes to reconstruct the top quark parameters (primarily, 
mt and Fj) from some exclusive jet observables, in presence of — > jet-jet 
decays. 

According to the HERWIG standard algorithm for PS, gluon radiation is 
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treated in the soft and coUinear approximation and no emission is permitted 

in the 'dead zones' (hard and large-angle parton radiation) seen in the left- 
hand side of Figs. 1-2. The current version of the HERWIG event generator 
has been improved by applying so-called ME corrections: i.e., the dead zone is 
populated by the use of the exact first-order matrix clement ('hard correction') 
and the 0{as) result is used in the already-filled region any time an emission 
is the 'hardest so far' ('soft correction')''. 




0.004 — 




log 73 



Figure 1: The t Wfb2g3 decay: (Left) The total (solid) and HERWIG (dashed) 
phase space in terms of the xi and x-j momentum fractions; (Right) The ys dis- 
tributions according to the new HERWIG implementation (dashed), to the old one 
(dotted) and to the exact 0{as) calculation (solid). Rates are at 360 GeV. 




,0 



x, log yg 



Figure 2: The e+e" -> tit2g3 production: (Left) Total (solid) and HERWIG (dotted) 
phase space limits for tt production in terms of the xi and x-^ momentum fractions; 
(Right) The ya distributions according to the old HERWIG implementation (dashed), 
to the new one, the latter populating either the small- and large-angle (solid) or only 
the large-angle region (dotted) of the dead zone. Rates are at 500 GeV. 
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A good measure to gauge the improvement due to the new algorithm can 

be gained from investigating the threshold value of the Durham jet-clustering 
algorithm ^ for all events to yield three jets, each with Et > 10 GeV and 
relative AR > 0.7 in transverse energy and cone separation. Just above ^/s = 
2m.i, where real radiation from the production stage is negligible, the beneficial 
effects of the new HERWIG implementation can be appreciated in the right- 
hand side of Fig. 1, in the case of top quark decays. At large ys values, where 
the fixed order 0{as) description^ is expected to well describe the underlying 
(hard gluon) dynamics, we see the latter to coincide with the current HERWIG 
description, which in turns retains sizable differences from the 0{as) result as 
2/3 (soft and coUinear gluon), where the all-order PS algorithm is known to 
be reliable. While similar improvements can be seen in the case of top quark 
production, it is in this instance more instructive to compare the same ys 
distribution as obtained from the old HERWIG description (incorporating ME 
corrections but not the complete ~ m^/s mass effects) to those given by two 
more advanced options, currently being implemented: the default one filling 
both the small- and large-angle dead zones and the intermediate one occupying 
only the large-angle region, both accounting for the above mass corrections (see 
the right-hand side of Fig. 2). The sizable difl[erences between the three clearly 
makes the point in favour of the more complete HERWIG implementation in 
view of future LC studies. 

Another example of ME corrections available in HERWIG is the case of 
multi-jet production in electron-positron annihilations via exact fixed-order 
e+e" — > n parton MEs, limitedly to the case n = 4. Here, the event generation 
can be initiated by the e+e~ qqgg and e+e" — > qqQQ MEs, with the four 
parton 'interfaced' to the subsequent PS (i.e., a hard cut-off is enforced in 
order to separate the partons thus preventing the MEs from diverging and the 
PS is attached to each of the latter only over the reduced phase space). This 
is an alternative to the standard description based on the e^e~ qq{g) hard 
scattering process followed by the traditional PS algorithm. While the new 
implementation is only reliable for n > 4 and for jet separations larger than 
the corresponding partonic cut-off, it remedies the drawback of the old one 
in describing the angular orientation of four-jet events, as seen by, e.g., the 
ALEPH Collaboration 11 (see also Ref. i^). 

The ability of a \IC event generator to describe correctly a four-jet final 
state is particularly relevant in the context of LC searches for a light Higgs 
state, predominantly produced via e+e^ ZH 4 jets events, as BR[Z 
jet-jet]=0.7 and BR[i7 bb] « 1. Take for example an hadronic sample 
composed of two &-jets and two untagged jets at a 500 GeV LC, selected by 
using the Durham jet-clustering algorithm with resolution 2/4 — 0.001 and with 
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soUd; HZ, dashes; 4j (HERWIG 5.9), dots: 4j (HERWIG 6,1) 
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Figure 3; Distributions for the polar angle of the reconstructed Z for the Higgs signal 
and the two description of the QCD four-jet background mentioned in the text. 

one pair of jets reproducing the Z mass within 10 GeV, i.e. |Mjj — Mz\ < 10 
GeV. Then compare the polar angle distribution of the di-jet system emulating 
a Z decay, for the QCD four-jet background and the above Higgs signal (for 
Mh = 150 GeV): see Fig. 3. It is clear that, in the region populated by the 
latter, the discrepancy between the old (HERWIG v5.9) and new (HERWIG 
v6.1) description of the former is rather large, similar in size to the signal 
excess. Needless to say, under these circumstances, a clear understanding of 
the behaviour of the QCD background is crucial in order to extract the Higgs 
resonance and measure its relevant parameters (Mh, Th, etc.) 

A more sophisticated implementation 'matching' MEs and PS, meaning the 
procedure of exploiting a combined approach which uses MEs for large phase 
space separations, PS in the infrared (i.e., soft and collinear) limit and modified 
(by Sudakov form factors) MEs in the overlapping region is available in 
APACIC-F-F, forn < 6 K 

2.3 Spin correlations and lepton beam polarisation 

Since v6.4, spin correlations are available in processes where SUSY particles, 
top quarks and r-leptons are produced Whenever these particles decay, 
2-, 3- and 4-body MEs are used to describe their dynamics (with or without 
the spin correlations). In the case of r's, an interface to TAUOLA is also 
available. Polarisation has been implemented for incoming leptonic beams in 

"Options similar to the one available in HERWIG for n = 4 can also be found in 
PYTHIA" and 4JPHACT". 

'A similar standalone implementation for n < 4 is available in HERWIG too 
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SUSY processes. These effects are included both in the production of SUSY 

particles and via the spin correlation algorithm in their decays. 

As an example of the impact of such spin effects, consider pair production 
of the two lightest neutralinos at an (un)polarised LC, via the following graphs: 




-X? 



I 

I 

^ > x^ 



with x5 the Lightest Supersymmetric Particle (LSP) and X2 decaying via: 




(a) Unpolarized (b) 6^6^ (b) e ^e^ 




1 2 3 1 2 3 1 S 3 

e^-p/rad e^-p/rad 

Figure 4: Angle between the lepton produced in e+e~ — > X2X1 ~^ ff XiXi if = e,^) 
and the incoming electron beam in the laboratory frame for ^/s = 500 GeV with: (a) 
no polarization, (b) negatively polarized electrons and positively polarized positrons 
and (c) positively polarized electrons and negatively polarized positrons. The solid 
line shows the default result from HERWIG which treats the production and decay 
as independent, the dashed line gives the full result from the four-body ME and the 
dotted line the result of the spin correlation algorithm. (The SUSY parameters are: 
Mo = 210 GeV, Aq = 0, tan/? = 10, Mi = 450 GeV, Afa.s = 350 GeV.) 

From Ref. we show Fig. 4, illustrating the angular behaviour of the 
above SUSY signals, as given by a typical MC implementation (HW), which 
does not include spin effects in production and decay and/or ME corrections 
in the latter, to the output of HERWIG v6.4 (HW+Spin). Once again, the 
improvement is clear if one confronts the latter with the distributions obtained 
from a complete calculation (ME) , wherein the event generation is based on the 
exact 2 — > 6 ME description. 

I thank the LCWS2002 organisers for the excellent atmosphere and stimu- 
lating environment that they have created during the workshop and The Royal 
Society of London, UK, for partial financial support in the form of a Conference 
Grant. 
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